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ABSTRACT 


The planning, layout and design of National Command 
Complex to resist all the nuclear weapon effects, to safeguard 
our political and military leaders in event of nuclear war has 
been illustrated 

An attempt has been made to combat the adverse effects 
of weapon by using Reinforced Cement Concrete as a basic 
construction material. The soil cover on top of the complex Is 
used to mitigate the radiation effects. 

An algorithm to design such a specialized 
specialised structure been suggested. A concealed and dispersed 
layout of the complex has been chosen to reduce the chances of 
success of the attacker’ s mission. 

The parametric studies have been performed to obtain the 
optimal depth of burial from blast, nuclear and thermal radiation 
considerations for various ranges of explosion. 

Static and dynamic analysis have been performed for onfe 
of the principal module of the complex using Structural Analysis 
Program - IV. The results are compared and a design force has 
been suggested. 

Finally, the module has been designed using Indian 
Standard Codes of Practice and checked for buckling stresses and 
critical moments 
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CHAPTER - I 


INTRODUCTION 


1 . 1 General 

Throughout the history man has sought to protect himself 
from hostile forces of his environment Apart from providing 
shelter from natural danger he has also striven to provide 
shelters from military threats. Seen historically, the need for 
defence and improvements in the means of defence were responses to 
threats and to counter-improvements in the method of attack 

1.2 Philosophy of Warfare 

From the earliest times, strongholds were solely 
defensives to resist attack by primitive weapons with no great 
range. The defenders needed an impregnable fortress, in which to 
shelter. The sophistication grew in response to increased numbers 
of deliveries from attackers, but the basic balance between 
attackers and defenders remained until it was upset by Romans' 
development of ‘siege warfare’ To counter this came construction 
of protective shelter inside the besieged area and the erection of 
tall shelters encircling walls. Thus the balance was maintained 
and see-saw battle continued until the end of nineteenth century. 

The growth in power and sophistication of explosives and 
fire arms changed the earlier philosophy during World War-I, in 
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favour of the attacker The defensive counter-measures merely 
increased the mass and strength of the structure in response to 
the more powerful attacks. The conflict until this time mainly 
confined to armed forces, who suffered twenty times more 
causalities than civilians The next major shift was due to 
advancement in aerial warfare widely used in World War II The 
civil populace behind the front line could for the first time be 
directly involved in the conflict. This led to further 
improvement in defence by development of public and private air 
raid shelters, which could withstand a local explosion and in some 
cases, even a direct hit 

1 . 3 Nuclear Environment 

The development of nuclear bomb and its first use at 
Hiroshima (15 kT) and three days later at Nagasaki (20 kT) 
radically altered the threat posed by the attacker The amount of 
nuclear weapon energy is so high that even if military targets are 
selected, the whimsically named collateral damage would result In 
an effect little different from that of such purely civilian 
targets as cities are chosen 

Parallel to this, came the remarkable improvement in 
accuracy of the delivery system of these devastating weapons, 
which is measured in terms of Circular Error Probabilities (CEPs). 
In fact some of the latest delivery systems have CEPs even less 


than a hundred meters. 
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1.4 Protective Shelters 


Thus in society an individual has only opportunity to 
safeguard an ordered civilization in long term, and that is by 
pressing the establishment of total disarmament The removal of 
fear of annihilation so that no nation has enough armament to 
attack any neighbor nation anywhere, must be the objective. Yet 
there are many powerful voices around the world, who still 
believe, ‘If you desire peace then be prepared for war’ This is 
dangerous nonsense when the weapons of war are not swords but 
nuclear warheads and it is a disastrous misconception to believe 
that by Increasing the total uncertainty one increases one’ s own 
certainty of survival 

Thus the only defence response made possible by these 
developments is to seek shelter in the immediate area to protect 
against the harmful effect of the weapons These shelters can 
broadly be divided into three categories from structural and 
utility point of view These are 

i) protective shelters for civilian population, 

ii) upgraded structures for industries 

iii) specialized structures (eg Hospitals, Water 

supply systems. Government buildings. Command 
complexes ) 

national security and chances of winning a 


As the 
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nuclear war are completely dependent on the continuity of the 
chain of the command at the highest level, this present study is 
focused on planning, design and construction of a Command Complex 
at the national level to safeguard our military and political 
leaders in the event of a nuclear war But the parametric studies 
and analysis of the present work are applicable for any hardened 
concrete structure to resist the nuclear weapon effects. 

1 . 5 Scope of the Present Study 

The scope of the present study is mainly to plan, 
analyze and design the 'National Command Complex’ A hardened 
concrete structure to mitigate the collateral nuclear weapon 
effects to the level tolerable to the housed 
operation/material/personnel is considered. To elaborate more, 
this protection entails 

i) nuclear blast and shock resistant design, 

ii) protection against radiation for specific 

period, 

iii) protection against thermal effects, and 

iv) shielding against other effects such as 

electromagnetic pulse (EMP) ,dust, debris. 


fires etc. 
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The aim of the present work is outlined as follows 

i) Construction of data charts for peak over 
pressure for various ranges, height of bursts 
and magnitude of explosion from available 
analytical expressions 

li) Establish pressure time histories for various 
possible cases of explosions, 
ill) Threat perception and Reliability analysis for 
dispersal of facilities. 

iv) Parametric study to obtain optimal depth of 
burial for various ranges, height of bursts 
and magnitude of explosion, taking soil 
arching effects into consideration 

v) Radiation analysis for various thickness, 
grade and type of concrete, 

vi) Study of attenuation of radiation with depth 
of burial for optimum depth taking bull tup 
factor into consideration. 

vii) Thermal Radiation Effect and its parametric 
study with depth of burial. 

viii) Static analysis of the principal shell module 
of concrete for most expected threat by using 
finite element method. 

ix) Dynamic analysis (Time History Analysis) of 
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the modules using Structural Analysis 
Programm-IV (SAP-IV). 

x) Design of the modules by Indian Standard 
Codes of practice 

1 7 General Standard Qualitative Requirements 

Some of the suggested qualitative requirements for such 
a complex are . 

i) Capable of housing approximately two hundred 

persons; 

li) Uninterrupted operation for a minimum of thirty 
days; 

lii) Alternate passages for facilities; 

iv) Sufficient dispersal and alternate emergency 

facilities; 

v) Concealed structure; 

vi) Emergency entry and exit doors, and 

vii) Internal and external communication 

1.8 Planning and Layout 

Based on the requirements listed in the foregoing and 
keeping in view the functional aspects, the proposed outline plan 
of the complex is as shown in Figure (Fig) 1.1. It caters for the 
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facilities as shown in Table 1 1 


Table 1.1 Suggested facilities of National Command Complex 


S.No 

Name of the facility 

Qty 

1 

Livlng/Off ices Rooms 

12 

2 

Administrative facility 

2 

3 

Power Plant 

2 

4 

Air Filiration Unit 

2 

5 

Hospital 

2 

6 

Cremation Place 

1 

7 

Fire Fighting Unit 

1 

8 

Water Tanks 

2 

9 

Decontamination Stations 

2 

10 

Security Rooms 

2 

11 

Chemical Toilets 

6 

12 

Central Hall 

1 

13 

Communication centers 

2 

14 

1 Meteorological Station 

1 

15 

Air Suction Unit 

2 


Apart from these the complex will also have passages 
to connect all the facilities for the movement of personnel and a 
Jeep size vehicle. 

Although the dimensions of these facilities may vary 
depending upon the actual functional requirements, for the purpose 
of this thesis, these have been kept at realistical values. 
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1.9 Organisation of the Thesis 

The present work has been organised into four chapters 
Chapter-I covers the planning and layout aspect of the complex 
keeping in view the qualitative and quantitative requirements. 

In Chapter-II the available literature on the subject 
matter has been reviewed in detail The aspects such as crater 
effects, ground motions and neutron radiation which are not 
considered in the present study have also been reviewed for the 
sake of completeness. 

Chapter-Ill deals with various parameteric studies 
performed from optimal dispersal, depth of burial and temperature 
considerations The latter part of the chapter contains analysis 
and design of one of the principal modules of the complex The 
results have been discussed at the end of the each section to 
facilitate better understanding. 

Chapter-IV concludes the present work with qualitative 
recommendations and future the scope of study for civil engineers 
in this challenging field . 



CHAPTER - II 


REVIEW OF LITERATURE 


2 . 1 Scope 

The literature review for the thesis has been classified 

as follows 

i) Blast and Shock effects of nuclear explosion 

ii) Crater effects, 

iii) Ground motions 

iv) Nuclear radiation effects 

v) Thermal radiation effects 

vi) Other effects (eg EMP,dust, debris and fires) 

vii) Threat perception and reliability analysis for 

optimal dispersal of facilities, 
vili) Soil arching relations and attenuation of over 
pressure with depth of burial 
ix) Finite Element modeling for shell structures 

x) Static and Dynamic Analysis 

XI ) Design of concrete shell structures 

2.2 General 

The nuclear weapon effects for protective construction 
vary with location and the magnitude of explosion. They can be 
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classified as air, surface or underground burst This location 
based classification is due to the variation in the quantum of the 
effects. A burst above 10m height is regarded as an air burst and 
the burst taking place below 5m is regarded as an underground 
burst [11 The relative effects due to various burst condition is 
shown in Table 2 2.1 [2]. It is evident that the near surface 
burst causes maximum damage in terms of blast and radiation It 
is more likely to be resorted by the attacker in the event of 
nuclear war. Hence, the analysis and design presented in next few 
chapteres caters for this worst possible condition. 

Table 2.2.1 Relative Degrees of Weapom Effects for Various 
Burst Conditions [2] 


Burst 

Condition 

Thermal 

Radiation 

Initial 

Hue 1 ear 

Radiation 

Ground 

Shock 

Air 

Blast 

Early 

Fallout 

Light 

Heal 

High Alt. 


** 

♦ 


* 


Air Burst 







Ground 







Water 







Confined 







subsurface 








2,3 Blast and Shock Effects 


When a nuclear weapon is detonated In the air, 
hydrodynamic expansion of bomb debris and radiation heated air 
produce a strong shock wave expanding in the surrounding air. Its 
interaction with ground surface affects the geometry and 
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characleristic of this expanding blast wave [1]. 

A surface or near-surface burst produces a hemispherical 
shock front that reflects from the ground The reflected signal, 
travelling faster than the incident shock in the heated air, 
eventually takes over the incident shock and forms a Mach stem 
that connects the ground with the triple intersection of incident, 
reflected and Mach stem shock to form, producing double peak 
phenomenon. This can be explained as shown in Fig. 2 3. 1 The 
incident wave (1) from a burst point above the ground is reflected 
in a regular reflection (2) until the angle it makes with the 
ground plane becomes large enough to allow the reflected wave to 
takeover and coalesce with the incident wave forming a Mach 
reflection (3). As the Mach stem grows, the contact surface (4), 
separating air that has experienced only the Mach shock from air 
that has been through both the Incident and reflected shocks 
trails behind and connects the triple point (5) with the surface. 
However, air flow above the contact surface does not completely 
turn parallel to the surface by the reflected shock and retains 
some downward momentum This results in eventual compression at 
the surface (6) as the downward moving double shocked air above 
the contact surface is stopped This compression (7) distorts and 
chews up the tail of the contact surface stream line (6) and rises 
the pressure. The Mach shock itself represents a lower pressure 
than that behind the reflected wave in the regular reflection 
region, but the second shock quickly builds to appreciably higher 
values at the transition. The important aspects of this 
phenomenon as shown in Fig. 2 3.2 (a) and (b) are as follows II] 




Fig 2*3,1 Kechanisni of Shock Wave for Near Surface Burst 



TJmt of!«r Oftowtiort 


Fig 2.3.2 (a) Variation of Over Pressure with 
Time at a given location 



Olttanc* from ConUr of Cipleilon 

Fig 2.3.2 (b) Variation of Over Pressure with 


distance at a given location 
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i) Time of arrival 

ii) Peak pressure 

iii) Positive phase duration 

iv) Pressure ^ Time History 

2.3 1 Over Pressure and Dynamic Pressure 

The transient pressure which is greater than the 
atmospheric pressure is defined as the over pressure. The peak 
over pressure is the maximum value of the over pressure at a given 
location, As the blast wave arrives at a given point, the over 
pressure rapidly increases from zero to the peak pressure. The 
magnitude of this peak over pressure and its variation depends on 
the location and the magnitude of the explosion. 

As the blast wave moves outward from the explosion, the 
flow of the mass of air behind the blast front produces a wind 
The resulting pressure is termed as dynamic pressure which is a 
function of density of air and the wind velocity behind the shock 
front The speed of the shock front U, is given by [1,3 and 4] 

r ^ P 1 

U = 1 + y - p— C2.3. 1 1) 

where , 

Assuming that air is an ideal gas 

y = Specific heat ratio =14 

= Ambient sound speed 

P 

o 


Ambient pressure 
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= Peak over pressure 

This approximation is correct within + 8 % for the peak 
over pressure less than 700 MPa [1] 

A finite time is elapsed between the detonation and the 
arrival of the blast front at a given location This time of 
arrival is dependent primarily upon the weapon yield and distance 
from the point of burst. 

The duration of the blast wave is characterized by two 
distinct phases, viz positive and negative phases. During the 
positive phase the over pressure rises very rapidly from ambient 
to peak value and then subsides more slowly to ambient pressure 
once again. For some height of burst in the mach reflection 
region, a second peak pressure increase may occur, which may be 
higher or lower than the first pulse As the distance from thp 
explosion increases the positive phase of the dynamic pressure 
becomes some what longer than that of peak over pressure. In the 
second phase a partial vacuum is created, then the wind blows 
towards the point of detonation. This phase is comparatively 
longer than positive phase and the peak rarely exceeds 28 
kPa. below the ambient value . Although smaller in magnitude it is 
important for design of blast doors and latches [1]. 

Yield Scaling 


The air blast parameters are normally given in terms of 
1 kt or 1 Mt. These parameters can be determined for other 
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yields of interest by scaling from the reference values [1] 


These scaling laws [1] are 

1/3 


f -(-S-) 

1 1 

LJ" fJL 

I w 
1 ^ 1 


1/3 


1/3 


(2.3 1 2) 


where 


R 

1 

W 

1 

R 

W 

t 

1 

t 

I 


Range at which pressure of interest occurs 

Reference yield 

Actual range 

Actual yield 

Reference time 

Actual timeand 

Reference impulse 


2.3.2 


Pressure - Time Histories 


Analytical expressions and several graphical methods are 
available [1 ,2 and 3] to find the over pressures and the time 
histones. These were developed and used since 1962 after the 
Nevada test sites results But the detailed experiments conducted 
In 1970* s have provided a better definition of the region of the 
double peaks, an important area in the early mach reflection 
region at high over pressures. These highly nonlinear analytical 
expressions to calculate the peak over pressures and corresponding 
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time history are given in Apendices ‘A’ and ‘B’ These 

expressions use the scaled ground range and height of burst and 
thus yield results corresponding to one kT Hence the different 
times and ranges must be scaled to the actual yield of 
requirements 

2.3.3 Approximate Analysis for Time Histories 

If computers are not available to solve these 

expressions then the time histories can be approximated by 
triangular equivalents, having the same initial peak-over 
pressures but different time depending upon the expected time of 
structural response. These equivalent triangular pulses are as 
shown in Fig 2 3 3.1 The durations are determined from the Fig 

2.3.3 2 . Reference [1] provides the procedural details for 

approximate analysis 

2.3.4 Dynamic Pressure 

In many a case the, drag forces associated with the 
strong transient wind behind the shock front are mainly 

responsible for actual damage, specially in the case of above 
ground structures. Some of the structures are described as drag 
sensitive, as they are quickly enveloped by the blast wave and 
thus the primary translatory forces are due to the results of the 
drag forces. These forces are the function of the size, shape and 
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transient dynamic pressure 

The analytical expressions for peak dynamic pressure 
dynamic pressure time histones, time of arrival, duration and 
impulse are given in References [1 and 31, As it is always 
advantageous to bury the structure for many obvious reasons hence 
the dynamic pressure is not of much importance in the present 
context However, this aspect is very important for the above 
ground structures and should be taken care while designing such 
structures The references [1,2] give details of this aspect for 
Lhe design of nuclear structures 

2,4 Crater Effects Due to Nuclear Blast 

This phenomenon occurs when there is very near, surface 
or under ground burst The approximate dimensions and other 
terminology associated with this aspect are as shown in Fig 2.4,1 

The dimensions are predicted from two tests conducted at 
Nevada test site, before nuclear test ban treaty, one in alluvial 
soil and other in hard rock. For other conditions these can be 
interpolated The scaling laws exist here to standardize for 1 kT 
explosion For the other yields crater dimensions can be 
predicted by applying a yield scaling exponent. 

The height of burst (HOB) or depth of burst (DOE) is 
measured from the original ground surface to the centre of weapon 
release. The most significant to catering is the low air burst, 
when HOB less than 3 m/kT^"^^. 


The volume of crater from air burst in various 
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homogeneous geologies are given in literature [1] and can be 
estimated approximately 

2.5 Ground Motions 

Another important phenomenon occuring due to blast in 
very near surface burst or ground burst is the induced ground 
motion as shown in Fig 2 5.1 Here the air blast loading is shown 
in super seismic condition i.e. the velocity of the air blast 
shock is greater than the propagation velocity of the peak 
compressive stress in upper material. The initial ground shock is 
a response to the air blast loading immediately above the point of 
interest. The ground shock component associated with the local 
air blast is said to be air slap induced. 

As the air blast shock front slows, the velocity becomes 
less than the propagation velocity of the peak compressive stress 
in the upper material The initial ground shock response at this 
time is caused by the air blast upstream from the point of 
interest Situations where the ground motions are outrunning the 
blast are of Importance from the analysis considerations [1]. 

The ground shock is divided in two components, the air 
slap induced and upstream induced components. The wave forms for 
each are superimposed temporally to arrive at the composite wave 
form. The ground shock from all sources other than local air blast 
is considered upstream ground induced shock, including shock from 
directly coupled energy and upstream air blast. This is further 
divided into body wave and surface wave components, because at 
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longer range the upstream induced motions are primarily associated 
with surface wave [1] 

To gain maximum advantage from the nuclear explosion on 
the strategical targets, the weapon is detonated in the air. 
Hence it very unlikely to expect any surface burst from the 
attacker Therefore this effect has not been considered for the 
design and analysis in the present context 

2 6 Radiation Effects 

One of the special features of the nuclear explosion is 
that it IS associated by emissions of the nuclear radiations, 
which is one of the governing criteria for design of such 
important structures [5] These radiations can be categorized as 


i) 

OL 

particles 

ii) 

P 

particles 

iii) 

r 

radiation 

iv) 

Neutron radiations 


2.6.1 (X and /3 Radiations 

The ranges and penetration power of these two radiations 
are very limited they cannot reach the surface of the earth from 
air burst even when the fireball touches the ground Hence the 
initial radiation consists of mainly y rays and neutron radiation 
within a minute of the explosion [2]. 
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Both of these radiations produce harmful effects on 
living organism as shown in Appendix (Apex) ‘C’ It is this 
highly injurious nature of these radiations combined with their 
long ranges that make them such a significant aspect of nuclear 
radiation. The energy associated with this is approximately three 
percent of the total energy release [1] 

2.6.2 y - Rays 

As it is evident that the y rays being the most 
important aspect of the nuclear radiation as for as the shielding 
IS concerned This is has been reviewed in greater detail under 
the following parameters. 

a) Sources 

In nuclear explosion the y rays are generated due to the 
fission process. They are also generated by the neutron produced 
in the fission process. When these neutrons are captured by non 
fissionable nuclei material they get converted into compound 
nucleus which is of very high energy state. The excessive energy 
may be liberated as y radiation called captured y radiation. Here 
the interaction of neutron with nitrogen nuclei is of Importance 
because y rays produced thereby is of very high energy and 
attenuates very little with distance. The y rays are also 
produced due to inelastic scattering As the fission fragments 
are highly radioactive in nature they also generate an appreciable 
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amount of the radiation in the very first minute of the explosion 
This radiation is known as delayed gamma radiation [2] 

b) Technical Aspect of y radiation 

Rigorous analysis of the radiation calculation due to 
rays needs the knowledge of interaction of the y rays with matter, 
as a result of which photons are scattered or absorbed These 
interactions can be classified as follows [5] 


i) 

Compton effect 

11) 

Photoelectric effect 

ill ) 

Pair Production 

Compton 

effect 


The total extent of compton scattering per atom of 
material with which the radiation interacts is proportional to the 
number of electron in the atom i.e to the atomic number for all 
materials, irrespective of their atomic weight. The effect is 
shown in Fig 2.6 2 l,(a] [5], 

Photoelectric effect 

It is an interactive process in which an orbital 
electron is ejected from an atom by a photon. This occurs when 
photon possesses a higher energy than binding energy of the 



25 


orbilal electron with which it collides All the energy of the 

incident photon in excess of the electron binding energy is 

transformed into kinetic energy of the ejected photo- electron 

As the binding energy of the orbital electron increases the 

photoelectric effects decrease rapidly as inversely proportional 

to the cube of energy, The photoelectric effect is proportional 

to the fifth power of the atomic number and the microscopic 

5 3 5 

cross section is proportional to (N Z E ) The other forms of 
radiation may also arise but they can be disregarded for practical 
purposes The phenomenology is sketched as shown in Fig 

2,6.2 1 (b) [5], 

Pair Production 


The process in which a photon is converted into a pair 
of electron, positron and negatron, in the coulomb field of an 
atomic nucleus is called pair production The incident photon is 
replaced by a posltron-negatron pair All the energy Is absorbed 
As the rest mass of an electron is equivlant to 0,51 Mev, Thus to 
have the kind of pair production incident photon should not have 
energy less then 1 02 Mev. The balance of the photon energy then 
appears as kinetic energy (KE) of the electron-positron pair. 

The probability of pair production is more when the 
atomic number Z of the absorber is high. The cross section can be 
approximated to Z . The cross section also increases with the 
photon energy, and at high energies it becomes the dominant 
interaction process of photons with matter (Fig, 2.6.2 1 (c)). 



Z6 
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(b) Photoelectric Effect 



(c) Pair Production 
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Hence it can be concluded that interaction of photon increases 
with the increase of atomic number of the material through which 
the y ray passes. Each interaction does completely absorbs or 
reduces the energy of incident photon Hence it can be said that 
material containing the elements of high atomic weight will be 
more effective as ^ ray shield Also with energy greater than 
1 02 Mev the absorption is minimum* This process of the 
interaction is shown in Fig 2 6,2 1 (c) [5 and 6] 

c) Photon Attenuation 

Consider a slab of absorber material placed between a 
narrow collimated source of mono magnetic photons and a narrow 
collimated detector as shown in Fig 2.6 2,2. The arrangements are 
such that the scattered radiation from striking the detector and 
any interaction of photons with the absorber material are 
prevented leading to complete removal of photon from the beam. 

Then, in a thickness dx of the absorber the fractional 
reduction in the photon intensity J is given as -dl/I, which is 
related to the narrow beam attenuation p and the thickness dx by 
the following relation 

dl/I = pdx (2,6 2,1) 

If the initial photon intensity is I^and the thickness 
of the absorber is x then 
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I M Jdx (2 6 2 2) 

I 

o 

and 

I = I e (26 2.3) 

o 

These equations are referred as exponential absorption 

equations 

The linear absorption coefficients are dependent on the 
density of the material The density often does not have a 

unique value but is dependent on the physical state of the 
material Thus to obviate the effects of variations the linear 
absorption coefficient for reference purposes is expressed a? 
mass attenuation coefficient (fi/p) [5 and 6], 

d) Photon Absorption 

The energy absorption coefficient are related with the 
loss or absorption of energy as a result of photon interaction 
with matter, whereas the attenuation coefficients concern the 
number of photons which are affected. Thus the amount of energy 
absorbed determines the effect of radiation Similar to the 
attenuation coefficients the values of the absorption coefficients 
are available in references [7J 
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e) Broad Beam Attenuation and Buildup Factors 


The above mentioned equations are only applicable to 
cases in which the scattered compton interaction may be regarded 
as essentially removed from the y ray beam, It reasonably holds 
good for narrow beams or for the shields of narrow thicknesses but 
fails for broad beams or thick shields, which is the case in 
nuclear explosion The allowance of multiple scattering is made 
by the provision of a 'buildup factor as shown in Fig 2 6.2.3. 
Hence these equations can be modified as [5] 

I = I B(x) e ^ 

O 

where 

B(x) Buildup factor 

The buildup factor depends on many factors such as 
energy of the photon source, angle of incidence, nature, thickness 
of the shield and the type of detector response The values of 
the buildup factors' for various ranges of energy and shield 
material is given in references [71. However, some of the relevant 
ones for the purposes for the thesis are given in Apex 


f ) Transmission of 'y ray and Source Geometry 

A point source is regarded as a volume of the source 
material the dimensions of which are small in comparison with the 
distance between the source and the position of the detector If 
a point source emits particles in all the directions 
isotropically, the radiation is assumed to be evenly distributed 
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2 

over the surface of a sphere whose area is 47 rx , where x is the 
distance between the source and detector The flux at a distance 
X IS given by [2 and 5] 

S 

(l> ~ — ^ (2 6 2.5) 

4 71 

where 

(j) = The particle flux 

= Strength of the source 

If it IS considered that the point source Is surrounded 
by a spherical shield then the flux density at a distance x from 
the source is 


= 


-tit 


4 tr X 


where 


( 2.6 2 6 ) 


t Thickness of the shield 

g) y Ray Spectrum 

The information concerning the energy spectrum of thp 
initial radiation is important because of the susceptibility of 
the living organism under various energy spectra. The attenuation 
properties of the air, shielding material and the response of the 
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detector instrument Is dependent on the energy distribution of the 
^--rays Although the exact spectrum at some distance from the 
source is difficult to obtain but with the knowledge of initial 
radiation spectrum at the source and existing data on scattering 
and absorption of different energy in the intervening medium can 
be employed to provide the spectrum at a given distance 
Reference [2] provides relevant information on this important 
aspect of the nuclear explosion. 

2.6.3 Neutron Radiation 

As in the case of the photons, which interacts with the 
electrons in atoms, the interaction of the neutrons with electrons 
IS negligible, They interact rather with the nuclei of the atoms, 
and their interaction is very much dependent on the KE of the 
neutrons. 

a) Nuclear Reactions 

Nuclear reaction generally be expressed by the 
following type of the equations [5] 

a + X ~> b + Y (2.6.3. 1) 

where 

X Target nucleus 


Y 


Resulting nucleus 
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a The bombarding particle 
b Emitted particle 

In the transmission of the neutron through matter many 
types of reactions or Interactions are possible In the shielding 
calculations for neutrons the reaction which are significant are 
[3] elastic, inelastic and neutron capture. 

b) Effective Removal Cross Section 

Most neutrons which penetrate thick shields are fast 
neutrons, having energy greater than 7 Mev. The effective removal 
cross sections are the cross sections which may be used in 
exponential absorption equation to obtain valid results concerning 
the attenuation of fast neutrons [5] Hence, 


I 


-I 


R 


( 2 . 6 . 3 . 2 ) 


Where 


^ Effective macroscopic removal cross 

section for material 
I Initial fast neutron intensity 


I 


Emergent fast neutron intensity 
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But, if the sufficient number of hydrogen atoms are not 
in the shield then the intensity of the emergent neutron may be 
higher than the calculated from equation 2.6 3 2 [5] 

Similar to the photons the energy spectrum for the 
neutrons from explosion of nuclear bombs are available in the 
literatures [2] The calculation for the integrated neutron flux 
can also be done in similar fashion 

For greater effects in terms of radiation special 
Neutron Bombs are available with most of the developed countries 
but they do not generate significant amount of shock effects 
The most part of the energy of these weapons are consumed to 
liberate radiation particles producing very little shock effects 
Since, the shock effects are the principal criteria for design of 
such structures, the effects of neutron radiation due to fission 
bombs in calculations can be neglected without much loss of 
accuracy IS] 

2.7 Thermal Radiation 

Thermal radiation is rarely a driver in the design of a 
hardened structure. By making provisions for protection against 
sterner threats of nuclear blast and nuclear radiations the 
designer almost automatically protects the structure against 
thermal radiation threat. However, it is pertinent to evaluate 
the thermal effects to the internal environment of the facilities. 
Also the study of thermal effects are important for the structural 
engineer to incorporate the temperature effects in the analysis 
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and for the selection of appropriate material accordingly. The 
approximate radiant exposures for ignition of various materials 
for low air bursts are given in literature [1 and 2]. 

2.7 1 Technical Aspect of Thermal Radiation 

Although spectroscopic studies made in course of weapon 
tests revealed that it does not behave exactly like a black body, 
the assumption of a black body behavior for the fireball, however 
serves as reasonable approximation in Interpreting the thermal 
radiation emission characteristics. For a black body the 
distribution of radiant energy over the spectrum can be related to 
the surface temperature by Plank's equation [2 and 8] as 


£ X 


8 n h 

c 


X 


5 


1 


h / \kT 
c 
e 


1 


(2.7.1) 


where 

£ Energy of the black body 

c Velocity of the light 

h Plank quantum action 

T Absolute temperature 

hc/X The energy of the the photon wavelength. 

The rate of emission of energy can be given for a 

A 


black body for a given wavelength[2 and 8] 
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= -^ £ (2.7 2) 


where 


J. Energy per unit area 


According to Stefan Boltzman law the total amount of 
energy J radiated per unit area per second is given by 


J = O' (2.7 3) 


where 


O' Stefan Boltzman constant 

= 2 15 

2 

= 5.67 E-5 erg/cm . s 

Next important parameter to be considered is the total 
thermal energy received at some distance D from the point of 
explosion which can be given by the law of physics [2] as 


Q = 


£ 

tot -kH 

.e 

4 Ti 


(2.7.4) 


where 
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k Absorption coefficient of the air 

E The fraction of the total weapon energy 

to t 

available as effective thermal radiation 
Q Total energy received at distance D 

it IS also important to know the total duration of the 
effective pulse for the calculation of the incident energy at any 
distance on the earth surface. Although no theoretical formulas 
are available for this time pulse time calculation, the tests 
conducted at the Nevada test site before nuclear test ban treaty 
provide some guidelines to evaluate this pulse time The graphs 
are available for various energies of the explosion in literature 
[ 21 . 

The attenuation of the thermal energy by a transient 
source in any medium Is a complex function of the conduction and 
radiation phenomenon It depends on the thermal properties of 
that medium and can be evaluated for the buried structure. 
Reference [9 and 10] provides the fundamental approach in this 
aspect of thermal radiation. 


2.8 Other Effects 

The other important effects of nuclear weapons which 
could be of importance to a design engineer are as follows 


i) Electromagnetic Pulse Effects 
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li) Debris 

iii) Dust and 

iv) Fires 

2.8.1 Electromagnetic Pulse Effects 

The nuclear explosion releases large amount of energy, 
Approximately one percent of it get converted in EMP. The nature 
of the pulse depends on the type, magnitude, range and the height 
of explosion. Its effect on any system depends on the physical 
and mechanical characteristics and the location of the system. 
These effects are more relevant to electronic equipments which are 
to be housed inside the complex Special care should be taken for 
installation of these equipments [1 and 2] 

2.8.2 Debris 

Debris is any material that can be carried along by the 
blast from a nuclear explosion. Such a material can cause 
causalities and damage to the facilities as a result of high 
energy impact or accumulation, or both This effect may cause 
severe damage to above ground and mounded structures . But due to 
many obvious reasons these command complexes are either shallow or 
deep buried, where the effects of debris are none Hence this 
effect has not been discussed in great detail. However, 
literature provides adequate information to combat this threat 


[ 11 . 
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2 8.3 Dust 

The dust is produced due to the following two reasons 

i) Propagation of the blast wave from the 

point of explosion 

11 ) Thermal energy produced due to the blast 

This phenomenon is again pertinent to above ground and 
mounded structures and therefore has not been reviewed keeping In 
view of the present requirements {1] 

2.8.4 Fires 

Although shallow buried structures have very little 
probability of catching fires, But civil engineers should avoid 
use of fire prone materials in construction of such complexes. 
References [1 and 2] give details of the threshold energy of 
various materials of civil engineer's concern [1 and 2]. 

2.9 Threat Perception and Reliability Analysis 

The threat perception and accordingly reliability 
analysis depends on various strategical and classified 
military informations If the information on the attacker's 

objective and defender’s actual requirements are available then 
this can be performed for actual scenario. In the present 
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context, due to absence of the actual information this has been 
performed only for dispersal of the facilities by assuming 
fictitious capability of the attacker However some of the 
guidelines available in Reference [1] can be summarised are as 
follows 


i) Develop demand scenario 

ii) Identify the failure modes and criteria 

lii) Calculate demand applied in each mode 

iv) Construct series and parallel network of 

failure mode probabilities 

v) Allocate goal to each failure mode 

vi) Calculate for each failure mode, capacity 

required to satisfy allocated share the goal 
vix) Identify worst case mean capacity for each 

failure mode and, 

viii) Write, for each failure mode, deterministic 

design specification for identified worst 
case mean capacity 

2.10 Blast Loading and Soil Arching Effects 

The choice of siting the complex is basically a 
compromise between the strategical and technical requirements 
The complex can be above ground, mounded or buried. The type and 
the quantum of the force subjected on the structure for each of 
these choices are totally different Therefore the approach for 
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analysis is different for each of these choices. The loading on 
the structures for below ground structures have been described in 
next few sections 

2.10.1 Air Blast Loading on Below Ground Structures 

The below ground structures can be categorized as 
shallow or deep buried depending on the depth of the burial (DOB). 
This division is relative and depends on the span of the structure 
The classification is as given below [1 and 11], 

i) Shallow buried ( 0. 2L <= DOB <= 1 5 L ) 

ii) Deep buried ( DOB > 1 5 L ) 

where 


L Span of the structure 

The design of deep buried structures is primarily 
governed by the soil overburden pressure. The blast has very 
little effect on the structure. These structures have other 
peculiar problem associated with them. They are generally not 
economical and hence not preferred unless warranted. Therefore as 
a part of literature review these structures have not been 
considered , 

The shallow buried structures experience partially the 
blast load due to soil arching and at the same time they are also 
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influenced by the soil overburden The total effects are 
different al different depths Hence for ease of the analysis 
these can further be classified as [1] 

i) Shallow Buried 

ii) Surface Flushand, 

ill) Mounded Structures 

2.10.2 Shallow Buried Structures 

The shallow buried structures have a DOB over the roof 
or crown, such that 0 2 L ^ DOB 2.0 L. where L is the clear 
span between the supporting walls of a roof slab or the horizontal 
diameter of an arch or cylindrical structure The structures 
located at 0 s DOB s 0.2 L are treated as surface flush The 
loads produced by the nuclear explosion on the structure buried 
in a medium with finite shear strength are Influenced not only by 
structure/medlum interaction as the loaded structure deforms but 
also by the effects of propagating transient stress wave [1] A 
typical shallow buried structure is as shown in Fig 2 10.2.1. 

2.10.3 Soil Arching 

It is defined as the ability of the soil to transfer the 
loads from one location to the another in response to a relative 
displacement between the locations Thus it is solely dependent 
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Fig 2. 10.2.. 1 


Shallow Buried Structure 
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on the properties of the soil This can effectively redistribute 
the loads on a roof slab away from the flexible center of span to 
hard support points such as columns and structural walls This 
has also been practically verified for buried arch structures The 
effect of soil arching is to force a high mode compression 
failure. For very hard structures the soil arch may not form 
because the shear capacity of the soil may have exceeded, and no 
redistribution would have occurred Various analytical 

expressions are available to incorporate these effects [1, 12, 13 
and 14] But the most generalized equation based on a 

differential approach that produces realistic results is given in 
Reference [1] and rewritten as 


Pq 

PsQ 


exp r- 2 Ko ( tan4> )( a + 1 )bL 


, . (2,10,3 1) 

where 

Pq Stress applied at the top of the structure , 

Pflo Ground surface over pressure 

Ko Coefficient of earth pressure at rest 

0 Angle of Internal friction for soil 

L Span length of the structure 

a Factor relating length of a structure to its 
span length 

b Factor relating depth of burial to its 
span length 

Fig 2.10 3.1 shows soil arching as a function of depth 


of burial. 
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2.10.4 Shallow Buried Rectangular Structures 


Predicting loads on structures is further complicated by 
the fact that stress wave propagating downward reflects at the 
structure surfaces which are in the path of propagation. This 
reflected wave is released by one of the three mechanisms When 
the structure begins to move as a rigid body with free field, by a 
rarefaction wave travelling inward from the lower stressed field 
surrounding the buried structure or by tensile wave reflecting 
back from the concrete-air interface at the bottom of the roof 
[1]. At very shallow depths it is very likely that a relief wave 
will be reflected back from the ground surface Assuming that the 


reflected 

stress IS relieved first 

by reflections 

from 

the 

concrete 

- air interface at the 

bottom of the 

roof 

the 

approximation for analysis can be done 

as follows [1] 




(rr(t) = cr/jf(t) ^2-t/tdj 

fort<= td 




crr(t) = crff(t) 

fort> td 




where 


(2. 10.4, 

1) 


(rr(t) The stress acting at the roof at any 
time t 

(rff(t) The incident stress wave 
td =■ 12 D/C 

C Compression wave velocity 
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D Thickness of the roof surface 

The time t can be taken as zero at the time incident 
wave strikes the roof surface Also for high yield weapons it is 
reasonable to assume the incident pressure as ground surface 
pressure but for low yield weapons less than 50kT it is 
recommended to take soil arching into consideration. Thus the 
equations can be rewritten as 


by 


(Tr ( t ) =(rf / ( t ) 2-t/tdj fort <= td 
crr(t)=C« <rff(t) for t> td 

where 


(2.10 4.2) 


Ca Soil arching ratio 


The stress acting at the base of the structure is given 


cPb (t) = P v(t) . . (2 10.4 3) 

where 

p Density of the soil 

Compression wave velocity of the soil 
v(t) Velocity of the structure. 

An acceptable approximation to the stresses acting at 


the base can be given by 
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(Tb (t) = crff (t) C2 10 4 4) 

The stresses acting at the side of the structure can be 

given by 


0's Ct) = o'ff (t) [2 10,4.5) 

where 

Coefficient of the lateral earth pressure 
at rest 

2.10.5 Shallow Buried Shells 

The analysis of a shell structure is complicated by its 
geometry [1,15]. The multidirectional effects become important at 
the instant the incident wave strikes the structure. The process 
of the Interaction can be described as shown ]n the Figs 2 10.5.1 
Ca) and (b). 

As the Incident wave strikes the upper crown of the 
structure, the stress at the crown is assumed to be double and due 
to this, structure begins to moving down as a rigid body. Within 
a very short period of time the entire half portion of the shell 
experiences stresses due to motion of the cylinder into the soil, 
while the incident wave has not progressed far beyond the crown 
As the incident wave propagates downwards, reflected wave 
continues to occur on the upper half of the cylinder because the 
peak reflected stress is a function of the angle between the 
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incident wave and tangent to the cylinder at the point where the 
incident wave contacts the cylinder, as well as function of the 
velocity of the structure at the time the incident wave arrives at 
the point The load modeling on the structure is shown in Fig 
2 10 5 2 til. 

By idealizing the cylindrical structure as octagon the 
load in the vertical region A can be given by 

F (t) = 0.8 O' (t) R L (2 10. 5. 1) 

a a 

where 

F (t) 
a 

R 

The load in the region B can be given by 

F (t) = 0 For 0<=t<03R/CL 

b 

F (t) = 1 2 (Tb (t) R L For t >= 0.3 R /C 

b ^ 

. . ,.(2.10 5.2) 

The total resisting incident load is taken as 

F (t) = 2 p C v(t) R L 
r c 


Total vertical load in the Region A 
Cylinder outside radius 
Stress in the region A 


. ... (2.10.5.3) 
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Hence ignoring shear the equation of motion can be given 

as [1] 


M v(l) = F (t) + F (t) - F (t) (2.10 5.4) 

a b r 

where 

v(t) Accelaration of the structure 

M Mass of the structure 

Another simplified method to analyze these buried shell 
structures by approxmating the blast force and lateral pressures 
Is shown in Fig 2.10.5.3 [ 1 ], is suitable for use with single 
degree of freedom system analysis which gives reasonable results. 


2 10.6 Mounded Structures 

The loads expected on these structures are intermediate 
between the above and underground cases. The inclined earth 
boundaries reduce the drag and reflected pressure to quite an 
extent but even then it is very significant up to a slope of 1:4 
The geometric complexity and non-uniform loading associated with 
mounded structures are not readily translatable into simple 
interaction model. However the bounding loads and motions of the 
mounded structures can be estimated by considering the above 
ground, buried or surface flush extreme. Since there are other 


(a) ARCH 


KPq 



(b) CYLINDER 

Fig 2.10.S.3 Simplified loading of Burled Arches 


Cylinders 



54 


effects of explosion which are more dominant on these structures, 
this choice is not suitable for the complex. 

2.10.7 Surface Flush Structures 

When the DOB is less than 0 2 L the structure is 
treated as surface flush due to Its different response than 
shallow buried structures The load acting on the roof of the 
structure can be analyzed as follows [ll 

(T (t) = crff(t) [2 - t / t’ 1 For t <= t’, 

r [ d J <1 

<r (t) = (rff(t) For t > t 

r d 

. (2.10 7.1) 

where 

O' (t) Stress acting on the roof at any 

r 

time t 

t* Z / C or 12 D/ C whichever is 

d s L 

smaller 

Z Depth of the structure roof <= 0 2L 

8 

As far as pressures at side and base of the structures 
are concerned they can be taken same as in case of the shallow 


buried ones. 
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2 11 Analysis and Design 


The analysis of shell structures can be performed either 
as thin or thick shells [16 and 17] The thin shell analysis does 
not take into account the stresses and moments In the direction 
along the thickness of the shell assuming it to be of a very small 
quantity The thin shell analysis approach can be adopted in the 
present context without much loss of accuracy in the results As 
the structure and loading are symmetric about one axis 

The membrane stresses and bending moments can be 
calculated by the following general shell equations [17] 




( 2.11 1 ) 


(2 11 2 ) 


Txy 




(2.11 3) 



(2.11.4) 


C2 11 5) 


(2 11 6 ) 


where 

F Stress function which gives the in plane 
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stress when the bending is also considered 
w Deflection along z- axis 

F 

D Flexural rigidity = and 

1 

V Poisson’ s ratio 

The analysis of the shells can be performed using SAP-IV 
[181. This current version of SAP-IV has facilities for static 
and dynamic analysis of linear structural systems This program 
has proven to be a very flexible and efficient tool This is 
coded in FORTRAN-IV and can easily be modified as per the 
requirements The fundamental approach for static and dynamic 
analysis using SAP-IV are given in next sections 

2.11.1 Static Analysis 

The static analysis performed by the SAP-IV can be 
summer ised as follows 

i) A static analysis involves the solution of the 
equilibrium equations [18] 

K u = R (2.11 1 1) 

which is followed by the calculation of element 
stresses 

ii) The load vector R is assembled at the same 
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time as the structure stiffness matrix are formed. The 

solution of the equation is obtained by the large 

capacity linear equation solver SESOL This subroutine 

uses the Gauss elimination on the positive definite 

symmetrical system of equations In the program, the 
T 

L DL decomposition of K is used and hence 

V = R (2 11 1 2) 

Y = D L u (2 11 1 3) 

where the solution for v in the equation is obtained by 
the reduction of the load vectors, the displacement 
vectors u are then calculated by back substitution [18], 

ill) After the nodal point displacement has been 

evaluated, sequentially the element stress displacement 
matrices are read and the element stresses are 
calculated, 

2.11.2 Dynamic Analysis 

The dynamic analysis is performed to obtain the 
following parameters 

-Free vibration frequencies 


-Mode shapes 
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-Displacements and rotations 
-Stress and bending moments 

The approach of the SAP-IV for the analysis is as 
described in the following 

i) The mode shapes and frequencies are obtained by 

solution of the generalized eigenvalue problem [ 18 ] 

K 0 = M (2 11 2 1) 


where 


w Free vibration frequency 

<p Mode shapes 

ii) In the dynamic response analysis the solution 
of the equation of the equation [ 18 ] 

MO’ + C u+ K u = R(t) C2 11.2 2) 


where 


R(t) Can be a vector, arbitrary time varying or 

effective loads resulting from ground 


motion. 
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ill) In the present context of the thesis R(t) is 
an arbitrary time varying load due to a nuclear blast at 
a given range and depth of burial Hence the Response 
History Analysis by direct integration is performed using 
SAP IV 

iv) The program uses Wilson-e method, which is 
unconditionally stable. 

v) The Rayleigh damping is assumed, which can 
easily be taken into account in the form of mass and 
stiffness proportion and is given by the equation 

C=aM+pK (2 10 2.3) 


where 

a and /? are mass and stiffness proportional damping 
coefficients. 

2.11.3 Design 

The design of structure is performed by using 
following Indian Standard Codes of Practices 

13:456(1978) Plain and Reinforced Concrete [19] 

IS-2210(1962) Shell and Folded Plates[20] 

Blast Resistant Design[21] 


13:4991(1968) 
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The check on the limiting stress is given in l22] for 
plates and shell elements 

2.12 Comments on Literature Review 

It is evident from the foregoing review that the 
literature provides sufficient guidelines and approach for 
designing these structures The approach may vary depending upon 
the requirements of the defender The following points are 
considered for design of complex 

i) It has been experienced that the above ground 
structures are economical only for peak over pressures 
less then 170 kPa. But for below ground structures the 
depth of the burial is a factor which depends on various 
economical and structural load parameters, which is one 
of the findings of the thesis 

ii) The choice of below ground structure is due to 
the following reasons 

It fulfills the military requirements 
It does not experience dynamic pressure 
which IS generated due to the following 
wind. 

The radiation is taken care by the soil 
thickness which is an economical and 
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effective radiation shield material. 

Thermal radiation does not have significant 
effects on the below ground structures. 

lii) Although the choice of burial reduces the 

impact of the nuclear forces but it presents a 
complicated loading due to the following 

Structures experience static soil overburden 
pressure in addition to the blast pressure 
The load prediction on these structures is 
complicated due to the soil arching effect 
The transient stress wave propagation in the 
medium with finite shear strength is a complex 
phenomenon 

The structure is subjected to active or 
passive earth pressure 

iv) The shells are preferred because of their 
shape which makes the magnitude of the net loading 
smaller than that developed on the box shaped ones. In 
addition they (domes and arches) have a great vertical 
and lateral strength characteristics and may be of the 
most economical form .However added cost is needed in 
terms of forms. 

v) The cylindrical and spherical shells and their 
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combinations are a better choice, as the loading is 
symmetric. These can be analyzed as thin shells using 
standard methods for membrane and bending components 

Vi) The finite element modeling of the modules can 
be done by assuming the modules on the rigid foundation. 
The other nodes have degrees of freedom in all the six 
directions The loading may be either static peak over 
pressure or the time histories depending upon the choice 
of the analysis The openings can be modeled by assuming 
different element type and material properties at those 
positions 

vii) The thin shell analysis is performed since 
the shells do not have significant displacements and 
rotations in the direction of the thickness 

viii) The membrane stresses and bending moments 
can be obtained using SAP-IV [18] with suitable 
modifications. 

ix) The design of the shell structures should be 
performed by using 13-456(1978) by limit state method for 
the principal modules of the complex 

x) The check for critical stress due to buckling 
should be performed as per 13-2210(1962) code for shell 
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structures 

xi) The Johnasen's failure theory for concrete 
shells/plates has been applied in the present context for 
checking principal moments in the structures The 
Reference 122] provides adequate details on this aspect 



CHAPTER - in 


DESIGN AND ANALYSIS 


3.1 Proposed Algorithm 

Taking into consideration the available technical 
information on the subject matter and the defence strategical 
requirements, an algorithm shown in Fig 3 1 1 is suggested, to design 
a specialized structure as command complex, subjected to multifacet 
forces This suggested algorithm incorporates the various possible 
effects of nuclear explosions at appropriate stages Thus a basic 
guideline to design such complexes of national interest is suggested. 

The planning and layout of the complex in the foregoing has 
been done as per this algorithm. 

3.2 Parametric Studies 

Various parametric studies which are general in nature for 
any type of the nuclear blast resistant structures may be carried out 
using the foregoing algorithm stated in Fig 3.1 1. 

3.2.1 Dispersal Analysis 

This analysis has been performed to obtain the optimal 
dimensions and maximum dispersal of the complex keeping in view the 
aticker’s capabilities. The assumptions made to perform this 
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analysis are given in succeeding paragraph 

Assumptions 

Assuming that the attacker has capabilities to 

i) launch only one bomb at the complex 

li) obtain the correct location if the target 

size = A , 
cl 

iii) correctly aim at the target if the target 

size = A 

Cl 

Based on the above mentioned assumptions the problem is 
formulated as given below 

Formulation for Dispersal Analysis 

If 

the max plan area of the facility = A^ and, 

the area of the complex = A 

s c 


then, 


the probability of obtaining correct location for a 
target size A^. 

pool(A^) = A^/ if A^ < A^^ 

- 1 If Aj. >=A^j. ..(3 2,1 1) 
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probability of correctly aiming at the target size 

pca(A ) = / K if A_ < A . 

f f Cl f ci 

= 1 If Aj. >= A^j (3 2 12) 

correctly aiming the target once located 
poci(A^) » pca(A^). .. . (3213) 

Circular probable error of the warhead 
delivery with respect to l<nown aim point 
Circular probable error in locating the 
attacker's aim point 

then 


Probability of 
pcal(A^) = 

Also, if 
CEP 

r 

CEP 

e 


the expected value of the warhead - strike distance R 
is (Fig 3 2.1 1) given by [1] 

R = 1.064 sqrt(CEP^^ + ' (3214) 


and if the 

probability of 

the warhead 

landing in the 

radius L 

r 




p('Wl) 

= L 

r 

if L < R 
r 


p(wl) 

= 1 

if L >= R . . 
r 

(3 2.1.5) 


where 

L is the radius defined by the attacker for 

r 


its successful mission 


Expected Value (R) and Coefficient of 
Variation (dR) of Weapon Strike Distance for 
Nearby Aim point 
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If, the probability of the explosion is p(e), then the 
probability of a successful mission 

ps(A^) =: pcal(A^) p(wl)._N_ p{e) (3 2 16) 

where 

N 

p(e) 

and the probability of failure of the attacker's mission 
can be given as 

pf (A^) 1 - ps(A^) (3 2. 1 7) 

provided there is only one facility existing in the 
radius L 

r 

It is evident from the foregoing formulation that the 
probability of the failure of the attacker’s mission is a function of 
the plan area of the facility. 

Hence, it is desirable to minimize the size of the 
facilities and disperse the same to an extent so as to reduce the 
probability of success of the attacker's mission A computer routine 
was written to solve the equation 3.2.1 6 to perform the parametric 
study to correlate the size of the facility and the probability of 
the failure of the attacker's mission. The study is preliminary in 


number of delivery 

probability of the explosion of the 
warhead once launched 
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nature as more classified military information which is secret in 
nature is required for a reliable situation However, it provides a 
basis to assess the reliability of such complexes 

Results and discussion 


The results plotted as shown in Fig 3 2 12 are based on 
following assumed capabilities of the attacker 


CEPr 

=100m 

CEPe 

=20m 

Lr 

=300m 

N 

=1 


It is evident from the results that to have very high 
chances of survival the facilities should be made smaller in the 
size For the facility size of 30 sq m the probability of the 
failure of the attacker’s mission is 98% which is quite reasonable 
The dispersal of the facilities should be more than the attacker’s 
defined radius for a successful mission Therefore the important 
facilities should have dispersal of approximately 300m and less 
important facilities may have dispersal of 150 - 200m to keep the 
compactness of the complex. 

3.2.2 Parametric studies for Optimal Depth of Burial 
(Lower Limit) 

As it IS evident that the blast pressure decreases 
exponentially as the DOB increases [1] » whereas the overburden 
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pressure increase with DOB can be taken as linear for shallow depths. 
Therefore these two phenomenon combined together provide an optimum 
DOB for minimum over pressure^^ 

In the present context the optimal depth of burial from 
minimum over pressure consideration for different ranges, height of 
burst and magnUudes of explosion, has been obtained by considering 
the soil arching effects and overburden pressure through following 
relations 


The overburden pressure P increases linearly with 

^ over 

depth by the following relation (Fig 3 2.2.1) 

P = r . dob , . (3.2.2. 1) 

over B 

where 

Z Density of the soil taken into 

S 

consideration 
DOB Depth of burial 

The attenuation of Peak Over Pressure with DOB can be 
given by the following soil - arching equation which is based 
on a differential approach 


P 


BO 


q 


exp 


-2 Ko (tan0)(a + llbL 


aL 


^ ■ 


(3.2 2.2) 


where 


P 

SO 


Peak over pressure and can be obtained by 
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Pressure (kPo) 



DOB (m) 


Fig 3.2. 2.1 


Variation of Pressure with Depth of Buriol 
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the analytical expression as given in 
Apex 'A' 

Pq The attenuated over pressure at the DOB 

under consideration. 

Coefficient of earth pressure at rest 

^ Angle of internal friction of the soil 

a Factor relating length of the structure 

to span 

b Factor relating DOB to the span 

L Length of the structure 

Hence, the total Over pressure at any depth is given as, 

(DOB) = P (DOB) + p (DOB) (3.2 2 3) 
total ovr q 


A computer program is written to solve the equation 
3.2 2.3 to obtain the optimal depth of burial (for minimum total 
over pressure) for various ranges and heights of burst for parametric 
study The following parameters have been assumed constant in 
the program The values of these parameters are given as 


3 

3' = 1,6 kg/m 

<t> = 20° 

K = 0,8 

0 

L = 6.0m 

a = 2/3 L 

b = 2/3 L 

Magnitude of Explosion = 20kT 


The results of the studies have been plotted as shown in 
Figs 3. 2. 2, 2 ,3.2 2 3 and 3.2.2 4 for different feasible cases 
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and are discussed in succeeding paragraph. 


Results and Discussion 


The study is of general nature and provides the lower limit 
or (maximum) depth of the burial for shallow buried structures from 
total over pressure considerations 

It IS evident from the graphical results that for near 
surface bursts the values of optimal pressures and corresponding 

t 

DOBs, for low ranges up to 500m are quite unstable These ranges 
should be avoided for construction. For HOBs up to 500m the optimal 
pressure and DOBs for 500-1000m range are more or less same and 
presents realistic values for design. Fig 3.2.2 4. can be taken as a 
design chart for obtaining correct DOB and corresponding pressure for 
precise calculations 

3.2.3 Parametric studies for Optimal Depth of Burial 
(Upper Limit) 

Another criteria for deciding the DOB for nuclear 
structures is radiation analysis. It is learnt from literature [1,3] 
that though the radiation is not a primary consideration for design 
but cannot be overlooked Therefore the upper limit of DOB is 
governed by radiation. In the present context upper limit has been 
obtained by keeping y - radiation limit as 0.3 rem for personnel 
inside the citadel. 

The principal approach adopted to obtain the minimum depth 
of burial for various ranges, HOB and Magnitude of explosion from 
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radiation consideration for different thickness is of the ordinary 
reinforced cement concrete is explained in Fig 3.2 3 1 The 


formulation of the problem is in 

the 

following 

Nuclear Fission Data 



1 

kT 


2 6 X lO^^Mev Energy 

1 

Fission 


200 Mev 

1 

Fission 


7 Mev release of energy 

1 

rad 

= 

0 01 J/Kg = 100 ergs/grm 

1 

Mev 

rr 

1 . 6 C-6ergs 


Assumptions 

i) Calculation of only primary and spontaneous 
flux IS considered It is presumed that if a 
person is able to sustain the initial 
concentrated dose, he/she would be able to 
survive under the reduced cumulative doses, 
due to residual radiation 
ii) The personnel living inside the citadel are 

exposed to the worst ^ ray flux, 
lii) For a common human being the exposure should 

be below 0,3 rem, 

iv) Residual radiation has not been considered as 

it IS an extremely small quantity 
v) The neutron radiation has been taken as half 
of the ^ ray equivalent, as the fission and 
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fusion bombs do not generate neutrons as 
primary weapon effects 

Formulation 


Let the flux at be A = d> 

Then 

$’ = C3 2 3.n 

A 

. 2 
4 TT y 

(3 2 3 2) 

B ■— — 

2 

4 71 AB 


“(ji / p ) BC 

= 4 e ^ ^ (3. 2. 3. 3) 

C B 

-(ji /p ) . t 

4 = 0 . e ® ° (3.2 3.4) 

D c 

where 

<l> y - ray flux at any point 

Mass absorption coefficient 
Pg Density of the soil 

Density of the concrete 
BC Direct thickness of the soil 

t Thickness of the concrete 


The total number of y - rays deposited in the human body 
inside the citadel can be expressed as (Worst possible 
case) 


y 


enh 


h A 


h 


(3 2.3.5) 
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where 


The average projected area of an average 
human being 


L 47ry^ 

(3 2.3,6) 

Hence, the energy deposited per unit weight of the average 
person (ergs/grm) 

r ,/Wt (^enh) .A * 1.6 E -6 .. (3,2 3 8) 

enh — Tf 

When quality factor (QF) = 1 then, 


-(M /p ) BC 


)] 


e^ 


y ^ (^enh) A • 1 6 E -6 . . (3 2 3 9) 

W(> 


If the buildup factor for soil and concrete are B^and 
B^respectlvely as calculated by Taylor's method for 8 Mev y ray 
energy (which could be the most expected average energy of the y 
ray emission in the nuclear explosion then, 


B = A 

a s 


-a fit 

.e ^ + (1 - A ) 

S 




(3 2.3.10) 


and 
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B=A e‘'+{l-A)e (3 2 3 11) 

c c c 

As, Ac, ots, ac, /3s, pc are Taylor* s coefficients 
obtained from the Apex ‘D’ 

Hence, the modified equations for flux after taking Buildup 
factor in to consideration are 

'^o » 2 3 12) 

(3.2313) 

A computer program has been written to solve above 
mentioned equations iteratively to obtain the minimum depth of burial 
for various ranges, height of bursts and thicknesses of the concrete 

The results plotted are as shown in Figs 3 23. 2, 3233, 
3 2.3 4 and 3 2.3 5. The results have been discussed in the 
succeeding paragraph 


Results and Discussio n 

This parametric study provides the upper limit of the depth 
of burial. It is necessary to bury the structure up to this minimum 
limit to avoid the ray exposure, which may cause the harmful 
effects to the occupants of the complex. 

This study revels that if DOB is more than 6m the structure 
is completely safe for near surface bursts. But DOBs of 2-6 m is 
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sufficient for most possible cases The concrete is no better than 
soil as for as radiation absorption is considered and therefore 
minimum concrete thickness from over pressure considerations should 
be taken for the analysis However to start the calculations for 
radiation a thickness of 250mm may be taken as standard 

3.2.4 Pareunetric study for Thermal Radiation 

Though the thermal radiation does not effect underground 
construction but it is important to ascertain the safe DOB from its 
effects The parametric study performed in the following provides 
useful results on this aspect The following assumptions have been 
made for the analysis 

Assumptions 

i) Nuclear explosion source is considered as a 

perfect black body, 

ii) No atmospheric attenuation is considered. 

Ill) Only one third of the total energy of the 

explosion is associated with direct heat 
radiation 

iv) The test results of the Nevada test site, 

conducted by the United States Department of 
Defence (USDD) has been taken to obtain the 
pulse time and other specific details for 
formulation of the heat radiation. 
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Formulation 


The fundamental approach adopted to calculate thermal 
radiation is as given in Fig 3241 

The approximate duration of the thermal pulse can be 
obtained in terms of tmax by the graphs as shown in the 
Fig 3,2 4 2 

where 


t = 0 032 W s (3.2 3 1) 

max 

W weapon yield in kT 

The total energy E emitted in the pulse can be 

tot 

approximated as 

Etot = Ji, (3 2.3 2) 

If there is no atmospheric attenuation the thermal 
energy Q on the earth surface at a distance D from the 
point of the explosion can be given by the following 
relation. (Fig 3 2 4.1) 

£ 

Q ~ tot -kD 
e 

4 TT 


C3.2.3 31 
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where 


k Absorption coefficient averaged over the whole 
spectrum of the wavelength 

A more useful formulation for thermal energy 
calculation can be given as 


fry 

Q = (3 2 3 4) 

4 n 

where 

f The fraction of the total yield in the form of 
the radiation 

T Transmxttancey a complex function of the 
visibility, absorption and the distance Its 
value can be obtained from Fig 3243 

The temperature at the earth surface due to the pulse 
source is governed by the following relation 
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I »-■=> [„ t] 
TT k ^ J 


1/2 


lerf c 


'^2(at} 


1/2] 


^ lerfc 


r (z\ A^) ] 

t 2 i 


where 

I 

A 

c 

k 

a 

t 


( 3235 ) 

Intensity of the incident thermal energy 

Area of the source 

Emissivity of the material 

The rma 1 conduc 1 1 v 1 1 y 

Thermal diffusivity of the material 

Pulse time 


The radiation of this transient temperature with time 
in atmosphere is given by Stefan - Boltzman law 


5Tg 

6T 


(P € 




(3 2 3 6) 


where 

<r Stefan-Boltzman Constant f 5.56 Watt/m ) 

e Emissivity of the soil 

Tg(t) Ground temperature in Kelvin at any time 

Te Environment temperature 

Z Depth under consideration 


The conduction at any depth Is given by the principle 
of heat conduction through infinite solid medium and 
can be expressed by the following relation 
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Tb - T(z,t) ( Z 'I C3 2 3 7J 

— = erfc 

Ts - Te >^2 

where 

Ts Temperature at the ground surface at any 

time t in Kelvin 

Te The environment temperature (K) 

T(z,t) Temperature at any depth z at any time t 
a Thermal diffusivity of the soil 

A computer program is written and used to obtain the 
temperature histones at various depths using above mentioned 
formulation for different ranges & HOBs, for a specific composition 
of the soil (Apex 'E*). The program assumes the values of some of the 
parameters as given below 

T = 0.9 

f = 0.33 

e = 0.93 

k = 0 13 W/m.deg 

W = 20kT 

The Results obtained by varying the other parameters in the 
equations 3.2 3 6 and 3. 2. 3. 7. The results obtained are plotted as 
shown in Figs 3.2.4. 4, 3. 2. 4. 5 and 3.2 4.6. These are discussed in 


following paragraph 



Maximum Temperoture (K) 



DOB (mm) 


Fig 3 2.4.6 Maximum Temperoture vs OOB 
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Results and Discussion 

The maximum temperature at any depth may cause thermal 
expansion of the structural material and is of interest from stress 
analysis point of view 

It IS evident that for near surface bursts the temperature 
on the surface may be as high as ZOOO^-'SOOO^c but it decays 
exponentially As the heat radiation from the top surface is faster 
than its conduction through soil layer Therefore the surface 
1000mm below ground virtually remain unaffected As other forms of 
radiation requires more DOB for the optimum choice » therefore the 
thermal radiation is rarely a driver in design As a rough 
approximation it can be stated that the DOB more than 1 m is 
sufficient to protect against the heat radiation for near surface 
bursts 


3.3 Pressure - Time History 

the pressure time history at the ground surface at any 
distance from the point of explosion can be obtained from analytical 
expression as given in Apex This provides a scaled output for 1 
kT of explosion Therefore the output available needs to be 
de-*scaled for the magnitudes of intrest. 

The arrival time and the duration of the positive phase can 
also be interpreted from the above expression. The value of peak 
over pressure as input for this expression is calculated from the 
expression given in Apex ‘A’ The time histones at various depths 
depend on the soil and structure Interaction with the blast wave As 
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a first approximation the time histones at any depth can be 
calculated with following assumptions 

Assumptions 

i) The load produced by ground shock on 

structures buried in a medium with a finite 
shear strength are influenced by structure 
medium interaction as the loaded structure 
deforms, besides by the effects of the 
propagating transient wave. 

ii) The reflected stress is relieved first by 
reflections from the concrete air interface at 
the bottom of the roof. 

Formulation 


With these assumptions the loading - time history at any 
depth can be formulated as follows 

The pressure time history at any given slant range 
from the point of explosion can be evaluated by the 
analytical expression as given in Apex 

For X >=xm and y <= 0.38 

P (x,y,t) =P (r, z) . (1+a) . (b. v+c) (3,3,1) 

s s 


and, for X < xm or y > 0,38 
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P„(x.y.t) = P (r.z) (b) (3 3 2) 

^ s 

where 

P^(x,y,t) Pressure at any time t for any 

range r from the point of explosion 
P^(r,z) Peak over pressure for any range r 

from the point of explosion and can 
be obtain by the analytical 
expression as given in Apex 'A* 
the other variables are as explained in Apex 

The stress at the roof surface at any depth without 
considering the soil arching relation can be obtained as 

cr t) = P (x,y,t) (2 - t/tj for t <= (3 3.3) 

r s d d 

O' (t) = P^Cx.y.t) for t < t . (3 3 4) 

r s Q. 

If the soil arching is considered then 

O' (t) = P (x,y,t). (2 - t/t^) for t <= t, 
r s a d 

<r (t) = C, P (x,y,t) for t < t (3.3.5) 

r a s a 

The stress acting on the side of the shallow buried 
structure can be given as 


(3 3.6) 
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If the soil - Arching is considered then 

<r (t) = C K P^(x.y.t) (3 3 7) 

a O S 


where 


1 ~ Sin <l> 
1 + Sin^i 


(3 3 8) 


This formulation provides vertical and horizontal pressure 
time histories acting on the roof and the side of the structure at 
the desired depth for shallow buried case As the overburden 
pressure on the structure constantly acts on the structure, the total 
vertical <r (t) and horizontal ir (t) histones can be 

rtotal stotal 

obtained as 


(3.3 9) 


(3.3 1C) 

where 

P = y . DOB 

vlro 


0- (t) 

r tota J 


O' (t) 
r 


vlro 


stotal 


k X (T ^ (t) 

0 rtotal 


A computer program has been written to obtain the loading 
on the roof and the side of the buried structure for various ranges 
from the point of the explosion The time-history obtained for a 
specific range (distance 300 m and HOB 100 m) and depths are shown in 
Figs 3.3.1. and 3.3.2 (a) and (b). These have been further used for 
dynamic analysis of the proposed modules In addition to the values 


Pr*i«uni {l(Po} 


: 

frO0C>0 “ 
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Rg 3 3 1 Surfoce Pressure - Time HjMory 


Pr«*9ure(kPfl) 


o 



Rg 3 3 2 (q) Verticol Pressure Time Hjslory at 
Vorious Depths 


Prp »8urp(kPo) 



Rg 3,3 2 (b) Horizonlol Pressure Time flt 

Vorious Depths 
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3f parameters assumed in previous sections the program assumes the 
i/alue of following parameters as 

t^ =12 D / C 

C = 600 m/s Compression wave velocity 

D = 250 mm Thickness of the roof 

3.4 Analysis 

The static and dynamic analysis of one of the the principal 
module of the complex has been performed 

The dimensions and other specifications of the module are 
as shown in Fig 3 4.1 

The analysis has been performed for 20 kT explosion at 300m 
away at 100m HOB 

M 20 concrete with density 24kN/m^ has been used for 
analysis purpose The other assumptions of parametric studies are 
valid for this analysis, 

3.4.1 Load Modeling on the Principal Module 

The load modeling of the principal module has been shown in 
Fig 3.4.2 (a) and (b) for static and dynamic analysis cases The 
enveloped view of the module has been shown with loading on the 
elements . 


3.4.2 Design Parameters 

Based and the parametric studies performed in the previous 
sections, the design parameters obtained for the static analysis 
which are given in Table 3. 4. 2.1 


100 


X 



Y 

Fig 3.4-1 Outline Diagram of the Principal Module 


Fig 3.4.2 (a) Load Modeling of the Principal Module for 


Static Analysis 



Fig 3.4.2 (b) Load Modeling of the Principal Module for 


Dynamic Analysis 
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Table 3.4.2. 1 Design Parameters For a Specific Case 


Range = 300m HOB = 100m Magnitude = 20 kT 


S.No 

DOB 

(m) 

Peak Over 

Pressure 

(kPa) 

Over- 

Burden 

Pressure 

9 -- Radiation 

(Rem) 

Thermal 

(K) 

1 

2 0 1 

352 3 

32 0 

0 3 

293 

2 

3 0 

280 0 

48 0 

<03 

293 

3 

4.0 

233 0 

64 0 

<03 

293 

4 

5.0 

202 

80 0 

< 0.3 

293 


Similarly the vertical and horizontal time histories 
acting on the structure at 300m range at 4m depth are shown in Fig 
3.3.2 (a) and (b). These histories have further been used for 
dynamic analysis 

3.4 3 Finite Element Modeling 

The suggested module of the complex has been modeled as 
thin shell elements of arbitrary geometry formed from four compatible 
triangles. The shell element uses the constant strain triangle. The 
central node is located at the average of the coordinates of the four 
nodes. The element has six degrees of freedom which are eliminated 
prior to the assembly, therefore the resulting quadrilateral element 
has four degrees of freedom i.e six degrees of freedom at each node 
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in the global coordinate system. 

The bottom nodes are constrained in all the direction and 
other nodes have degrees of freedom in all the six directions 

The openings In the module have been taken care of by 
providing different element type at the opening with negligible 
elastic modulus and the same thickness value to maintain the 
continuity and compatibility 

The thermal radiation effects on the structure due to blast 
have been taken care of by providing the temperatures at nodes 
according to the results obtained by parametric studies of section 
2 2 4 

The Global and Elemental coordinate system are as shown in 
the Fig 3 4 3.1. 

The finite element modeling of the module of the complex 
are as shown in Fig 3 4 3.2 

A preprocessor was developed to provide input to SAP-IV for 
static analysis The loads on the shell elements act normal to the 
surface for analysis. The stresses and deflections of all the 
elements are recovered in output which is of interest for the design 
The following parametric studies are performed to obtain the desired 
mesh size for optimal results 

3,4.3. 1 Parametric Studies 

The following parametric studies have been performed 

i) Stress on any arbitrary chosen element vs size 

of the mesh (Fig 3433) 

ii) Bending moment per unit run on any arbitrary 
element vs size of the mesh (Fig 3 4.3.4) 

iii) Maximum stress vs size of the mesh 
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(Fig 3 4 3 5) 

IV ) Maximum bending moment vs size of the mesh 

(Fig 3436) 

v) Processor time vs number of elements 

(Fig 3437) 

Results and Discussion 


The parametric studies performed above for optimal size of 
the mesh reveals that no of elements should be taken as 1600 - 2000 
(250x250mTti •~200x200mm) for optimum results if the thickness of the 
model IS uniform as 250mm, However, with increasing thicknesses at 
critical sections for stress reduction demands a fresh parametric 
study for the optimum results The element should be kept as square 
element. 

3.4.4 Static Analysis 

The parametric studies performed in the foregoing revealed 
that the openings and the joints of the shells are critical due to 
stress concentration Even for the vertical pressure of 200 kPa 
the stresses on these sections exceeded the critical stresses almost 
by 200%. Therefore the mesh size was readjusted with increasing 
thicknesses at these sections (ring beam at the joint of the shells) 

The parametric studies were once again performed to obtain 
the desired results It was found that the mesh size of 410ininx400miin 
(585 elements) provides the optimum stresses. 

The static analysis performed is performed with this 
readjusted mesh for various load conditions and varying thicknesses 
The results obtained are as shown in Table 3 4 4 1, (a) and (b). 
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Table 3. 4. 4.1 (a) Results of Static Analysis 


s 

DOB 

(m) 

Peak 

over 

pre 

saure 

(kPa) 

Thicknesa 

(mm) 

Max Stresses 

(kPa) 

Cyl Inder 

Dome 

Near 

opening 

0 

(20 ) 

Rest 

cyl 

Ring 

beam 

Rest 

dome 

S 

X 

5 

y 

S 

xy 

1 

2 0 

352.3 

600 

300 

400 

200 

0. 15E4 

0 36E4 

0 63E3 

2 

2 0 

352 3 

500 

250 

300 

150 

0 20E4 

0 41E4 

0 67E3 

3 

2 0 

352 3 

400 

200 

200 

100 

0.27E4 

0.52E4 

0, aOE4 

4 

3 0 

280,0 

600 

300 

400 

200 

0. 13E4 

0. 32E4 

0,62E3 

5 

3,0 

280.0 

500 

250 

300 

150 

0.16E4 

0 36E4 

0 64E3 

6 

3.0 

280 0 

400 

200 

200 

100 

0.21E4 

0 45E4 

0 7gE3 

7 

4 0 

233 0 

600 

300 

400 

200 

0 13E4 

0 30E4 

0 63E3 

8 

4.0 

233 0 

500 

250 

300 

150 

0.13E4 

0 33E4 

0 64E3 

9 

4.0 

233.0 

400 

200 

200 

100 

0 18E4 

0.40E4 

0.67E3 

10 

5 0 

202.0 

600 

300 

400 

200 

0 13E4 

0, 29E4 

0 63E3 

11 

5.0 

202,0 

500 

250 

300 

150 

0 13E4 

0 31E4 

0 64E3 

12 

S.O 

202 0 

400 

200 

200 

100 

0.15E4 

0.37E4 

0 67E3 


Results and Discussion 

The results obtained for various thicknesses of the module 
supports the fact that the critical sections are the openings (20° 
from the opening) and the joint of the shells. Therefore additional 
thickness are required at the Joint in form of ring beam and at the 
openings The compressive stresses are the governing stresses for 
the design of the modules. The maximum deflections in principle 
directions are as shown in Table 3.4.4 1 (b). 
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Table 3.4.4.1*(b) Results of Static Analysis 


s 

DOB 

(m) 


Thickness 

(mm ) 

Hay. Deflections 

(mm ) 




Near 

opening 

0 

(20 ) 

Rest 

cyl 

Ring 

beam 

Rest 

dome 

X 

transl- 

ation 

Y 

transl- 
at 1 on 

Z 

transl- 

ation 

1 

2.0 

352 3 

600 

300 

400 

200 

0 22E;-2 

0 24E-2 

0 22E-2 

2 

2 0 

352. 3 

500 

250 

300 

150 

0 31E-2 

0. 30E-2 

0 31E-2 

3 

2,0 

352 3 

400 

200 

200 

100 

0 52E-2 

0 59E-2 

0 50E-2 

4 

3.0 

280 0 

600 

300 

400 

200 

0 20E-2 

0 21E-2 

0 20E-2 

5 

3.0 

280 0 

500 

250 

300 

ISO 

0.27n-2 

0 30E-2 

0.27E-2 

6 

3 0 

280.0 

400 

200 

200 

100 

0.45E-2 

0 52E-2 

0 43E-2 

7 

4.0 

233,0 

600 

300 

400 

200 

0 18E-2 

0 19E-2 

0 18E-2 


m 

233.0 

500 

250 

300 

150 

0 25E-2 

0 27E-2 

0 24E-2 



233 0 

400 

200 

200 

100 

0.40E-2 

0.46E-2 

0.39E-2 

10 

5.0 

202.0 

600 

300 

400 

200 

0,17E-2 

0 18E-2 

0 17E-2 

11 

5,0 

202.0 

500 

250 

300 

150 

0 23E-2 

0 26E-2 

0.23E-2 

12 

5.0 

202 0 

400 

200 

200 

100 

0.37E-2 

0.43E-2 

0, 36E-2 


3.4.5 Dynamic Analysis 


As a first step to dynamic analysis natural frequencies and 
the mode shapes are calculated for three different specifications of I 

the structure. 

, I 

I 

A preprocessor is developed to provide input to SAP-'IV, I 

The first and second frequencies obtained for various cases are as I 

i 


shown in Table 3.4.5. 1. 
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Time History Analysis 


For dynamic analysis the same module of the comple>c is 

loaded with two time histones acting simultaneously These histories 
are the same which were obtained in section 3,3 One of these acts 
vertically (X-direction) on spherical shell (top dome) for 90s The 
other acts horizontally (normal to the surface) on the cylindrical 

portion of the module for the same time The horizontal time history 
has further been resolved in global Y and Z directions for the 

loading on each element. The force on each element has been provided 
with the help of time multiplier function in SAP-IV 

The time step has been chosen as 1 s. The damping has been 
assumed as 2%, 

A preprocessor has been developed to provide above 

mentioned details as formatted input to SAP-IV for time integration. 
The load modeling is as shown in Fig 342 (b) 

The maximum stresses, deflections and the corresponding 
time are recovered as an output are given in Table 3.4.5 2, (a) and 
(b). 


Table 3.4.5. 1 Natural frequencies of the module for 

various stiffnesses 


s 

no 

Thickness 

(mm) 

Frequencl es 

(Hz ) 


Near 

Rest 

Ring 

Rest 

First 

Second 


open 

cyll 

beam 

cyll 




0 

-nder 


-nder 




(20 ) 






X 

600 

400 

400 

200 

0. 3644E2 

0.4361E2 

2 

500 

300 

300 

150 

0. 3424E2 

0 3842E2 

3 

400 

200 

200 

100 

0. 2811E2 

0 3099E2 



112 


Table 3.4.5 2 (a) Results of Dynamic Analysis 


s 

DOB 

(m) 

Peak 

over 

pre 

ssure 

<kPa) 

Thickness 

(mm) 

Max Stresses 

(kPa) 

Cyl Inder 

Dome 

Near 

opening 

, 0 
(20 ) 

Rest 

cyl 

Ring 

beam 

Rest 

dome 

s 

X 

s 

y 

S 

xy 

1 

3.0 

280 0 

600 

300 

400 

200 

951.6 

1520. 0 

471 0 

2 

3 0 

280 0 

500 

250 

300 

150 

1189 0 

2030 0 

666 3 

3 1 

3.0 

280 0 

400 

200 

200 

100 

1641.0 

3015 0 

1094 


Table 3. 4. 5. 2 (b) Results of Dynamic Analysis 


5 

DOB 

(m) 

Peak 

over 

pre 

ssure 

(kPa) 

Thickness 

(mm) 

Max Deflections 

(mm ) 

Cyl inder 

Dome 

Near 

opening 

, 0 
(20 ) 

Rest 

cyl 

Ring 

beam 

Rest 

dome 

X 

transl- 

ation 

(m) 

Y 

transl- 

ation 

(m) 

Z 

transl 
ail on 

(m) 

1 

2.0 

280 0 

600 

300 

400 

200 

0.00142 

0 00113 

00134 

2 

2 0 

280 0 

500 

250 

300 

150 

0 00208 

0 00189 

00226 

3 

2,0 

280.0 

400 

200 

200 

100 

0.00370 

0-00203 

.00483 


R 

X 


S 

xstatlc ^ 1 T 

= g - 1 J 

xdyjiami c 


7 ®ydy...lo 


xystat Ic 

3 

xydynamic 


S 1 0 


where R , R R are the ratios of static and dynamic stresses 

X y xy 


Results and Discussion 

The comparison of the results of static and dynamic analysis 
reveal the fact that the static analysis for blast load is over 
conservative. Therefore the dynamic analysis of these structures are 
recommended for economy in design of the complex As the 
loading on these structures are one time measure, therefore 
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these structures should essentially remain safe from static 
overburden pressure 

3.5 Design of Concrete Shells 

Keeping in view the results of the analysis in the 
foregoing the design of the principal modules have been performed as 
per Indian Standard Codes of Practice [19,20 and 21] 

3.5.1 Design 

The moments of resistance of rectangular section based on 
the assumptions of 37-1 of 15*456(1962) for section without and with 
compression reinforcement Is given in Apex 'E" of the code has been 
used as a basis of the design The elements have been assumed as 
rectangular sections 

The element are checked for critical buckling stresses as 
per 15.2210(1962) The procedure is rewritten in the succeeding 
paragraph. 


3.5.2 Critical Buckling Stress 

Instability in a shell may be caused by 

- Local buckling in zones submitted to compressive 
stresses 

- Brazier effect 

^Combined effects of bending and torsion 
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Cylindrical Shell s 


The permissible buckling stress in a cylindrical shell may 
be calculated as follows 


0.25 F 

f = (3 5 2.1) 

perm . t 


F 

cr 

where 


F Cube Strength of the concrete at 28 days 

e 

F Critical stress determined in accordance 

cr 

with IS 2210(1962) 


Buckling in Doubly Curved Shells 

For spherical shells the permissible buckling load per unit 
area from consideration of elastic stability is given by the 
following relation 

P = 0.1 (3.5 2.2) 

perm ^ 2 

where 

E Elastic modulus of concrete 

d Thickness of shell and 


R 


Radius of Curvature 
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3.5.3 Yield Criteria 

The maximum moment at any corner of any element has been 
checked by assuming the elements as plate elements Johnasen’ s plate 
theory is thus applied to check the ultimate moments in each element 
The principle approach adopted is given below 

Let the reinforcement be placed in two orthogonal 
directions, one along radial X-axis and the other along the 
circumferential Y-axis of the shell. Since the reinforcement may not 
be same in both the directions, the corresponding plastic moments 
will also differ and the element will be orthotropic with regard to 
its ultimate strength. 

Then the ultimate bending moments can be denoted by M 

pX 

and N for the positive bending and M' and N'py for negative 
py ^ pX 

bending. 

Yielding occurs when the bending moment on the cross 
section with angle a reaches a certain value. Refer Fig 3.5 3 1 

Then according to Johansen’s failure theory expressions for 
lower Yj^ and upper yield Y^ curves can be given as 

Y MpOi = M , cos^oi + H .sin^a (3.5 3.1) 

L pX py 

Y,, Wpa = M' cos^a + W’ sin^a (3.5.3 2) 

U pX py 

The yield condition in terms of the components W^, and 

M IS given 
xy 


M = M cos^oi+ M sin^a + M.sin 2 a 
ax y xy 


(3 5 3.3) 
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Then the condition to be checked is whether 

M > M (3 5 3 4) 

pa a 

A common computer program has been written to check the 
buckling stresses and permissible moments at all the corners of each 
element applying the above principle As the buckling stress are the 
governing stresses therefore minimum reinforcement (ISmm 0 250 mm 
c/c) has been provided in both the principle directions The 
stresses obtained from the analysis are well within the permissible 
limits. 
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Fig 3. 5. 3.1 


Yield Criteria for shell Eleroenls 





CHAPTER - IV 


CONCLUSIONS 


4.1 Suggested Algorithm 

The algorithm suggested in section 3.1 incorporates all 
the effects of nuclear weapon which are governing forces for the 
design of the complex The planning, layout and design of the 
complex has been done as per this algorithm. 

4 . 2 Parame tri c Stud i es 

The parametric studies performed drive home the 
following points 

i) The probability of the failure of the 

attacker mission is a function of the dispersal of the 
facilities and the size of the modules Therefore the 


complex should 

be 

dispersed 

to the 

maximum 

and the 

dimensions of 

the 

facilities 

should 

be as 

small as 


possible 

li) Parametric study performed for optimal depth 

of burial from blast pressure considerations reveals that 
burying structure deep in the ground is not only 
uneconomical but it has to designed for higher soil 


overburden pressure 
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4.3 


The shallow buried structures are advantageous 
from blast, nuclear and thermal radiation considerations 
compare to above, mounded , surface flushed and deep 
buried structures . 

iv) The soil cover on the top of the structure Is 
is more economical than the thickness of the concrete 
from radiation point of view, 

v) Upper limit of burial should be based on the 
radiation considerations* 

VI ) Although the specialized concretes are 
available to combat the radiation but the ordinary RCC 
serves the purpose in present context 

vii) It is evident that the thermal radiation Is 
rarely a driver in design but the upper layer of the 
soil up to one meter depth is susceptible to the thermal 
radiation Therefore the structures should not be 
mounded or surface flushed 

Static and Dynamic Analysis 

i) The static analysis carried out on the 
principal shell module supports the fact that the 
openings and the Joints of the shells are critical from 
the stress concentration point of view* These should be 
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designed very carefully 

ii) The size of the mesh and its refinement with 

varying stiffness of the elements is important for 
^ desired results 

iii) The comparison of the static and dynamic 

stresses revals the fact that the static analysis for 
blast force as over conservative Dynamic analysis should 
be resorted with taking soil arching into consideration 
As the analytical expression for the pressure-time history 
is approximate, the dynamic analysis with a factor of 1.5 

IS recommended. 

4.4 Scope for future study 

The design of structures to resist against nuclear 
weapon effects is a challenging field for civil engineers The 
nations around the world have taken this aspect very seriously in 
recent years due to nuclear threats. 

We may have to construct new structures or upgrade 
existing construction to cater for these threats which 
automatically opens new fields for research The nonlinearity may 
be increased to obtain more accurate results for these effects, 
Also, the soil arching relations can be experimentally be found 


for realistic results. 


CALCULATION OF AIR BLAST PEAK OVER PRESSURE 


Given 

scaled range in km/ki 

1/3 

scaled Height of burst in jcm/kT 
Magnitude of explosion in ki 

Then 


X 

y 

m 


r 



2 


y 


) 


1/2 


2 = y/x 


Peak over pressure in kPa is given by 


Ps(r,z) = p 


■ 10 47 

- (a rf 


+ b(z) 

(a rf‘” 


. d(z) e(z 
+ 

1 + f(z) 


in which 


h(z,r,y) 


-(64.67 z^+ 2905) 

1 + 441.4 


1 389 z 
1 + 49.03 2^ 


2 

.0014(tf r) 

(1 - .158(a.r) + 0486 (a,r)^^ + 


where 


a = (0.3048)"\ kft(km) 


+ h(z,r,y) 

(« A"’ 


8 808 

+ 

1 + 154 5 z^'^ 


00128(a.rf ).(1 + 2 y) 
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P = (100/14. S04) , psi (kPa) 


a(z) = 1.22 - 


3.908‘'z 


1 + 810 2 z 


b(z} = 2.321 + ^ 


18 


1 + 1 113 z 


18 


03831 2 


17 


6692 


1 + 02415 z 


17 


1 + 4164 I 


c(z) = 4 153 - 


1 149 z 


18 


1 + 1 . 641 z 


18 


1.1 


1 + 2 771 z 


2.5 


d(z) = -4. 166 + 


25 76 z 


1.75 


1 + 1.382 z 


18 


8 257 z 


1 + 3 129 z 


e(z) = 1 - 


.004642 z 


18 


1 + .003886 z 


18 


f(z) = .6096 + 


2.879 z 


9 2S 


1 + 2 359 z 


14.5 


17 25 z 
1 + 71.66 z 


g(z) = 1 83 + 


5.361 z 
1 + .3139 z* 


' .I''--- 
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Apex ‘B* 


ANALYTICAL EXPRESSION FOR PRESSURE - TIME HISTORY 


For X >= xm 


and 


y <= 0 38 


P (x.y.t) 

9 


P (r, 2 ) 

B 


( 1 + a ) ( b V + c ) 


( 100 ] 
[ 14.504 J 


kPa 


For X < xm or y > 0 38 


Pjx,y,t) =p^(r,z) .( b) I 50 ^ ] kPa 
Where 

1/3 

X scaled ground range, kft/kT 

y scaled height of burst, kft/kT* 

1/3 

t scaled time msec/kT 


a = 



20 N 

20 

1 + e ^ 


d = 0 . 27 + 


58300y 


+ O 27 (e)- 0 . s - 


26667 + looooooy 


0 5830QOy 
26667+lOOOOOOy^ 


I X ^ xm I 
xe - xmT 


1 + 337 y 


+ 0.91 y 


3 039 
xe = 


1 +. 6.7 y 
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= u for X < xm 

j. xni xm 

ta - u + w ll ' — — for X >= xm 


2 3 8 

(0 543 -21. 8r + 386r ♦ 23B3r )r 


„ -10 2 -64 -56 28 

2 99 10 -1.91 10 r +1 032 10 r -4.43 10 r *(1 632+2 629+2 69r )r 


3.0137 


2 3 8 

(l 086 - 34 605r + 486 3r + 2383r )r 


-13 ,..-92 -6 4 -5 6, 2. 8 

10 -1.2128.10 r +4.128 10 r -1.116 10 r +(1 632+2 629r+2,69r )r 


dp 


(- 


1640700 + 24629 ta + 416 15 ta 1 
10,880 + 619 76 ta + 


rfo.4 . .-o-°qi 204 i iv^ . 

1 + 0.001559 


0426 + 


0 5486 ta'^^ 

1 + 0 000357 ti*® 


]-] 


= 1 _ 

1 + ].i.io^°y’^ 


r 2.441 10" V I ( 1 ] 

I 1 -+ 9 lo'V J ^ 4 41.10"“+ J 


:■= f 0 01 
U + 0 


01477^ ^ 7 402 10- = t.";° . 0 |_ , 


005836ta 1 + 1. 429,10"®tl' 


0.7076 - 


3,077. 10"^ta^ 

1 + 4,367 lO'^tl 


g = 10.0 + 


77.58 - 


64.99ta 


0. 12S 


1 + 0 04348ta ' ® 


)• 


(s) 
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h = 2 753 + 


0, 05601 ta 


1+1.473 


( — 
'■1+3 


01769ta 


207 10^° ta 


0 03209ta 


1 * 25 


1 + 9*914. 10 


-8 




(s) 


V = 1 + 


/ 116 s^ 

I 3 28 IQ y I [ ga j 

i . 12 6 75 I 3 ] 

'^1 + 15.10 y ^ ^ 6 13 + ga ^ 


1 


1 + 9 23 


2 

e 


f 240. 9x^ ] , .7-1 

1.04 - (ga) 

M+ 231. 7x ' 


^ 9 -» 

2, 3, 10 y 

(1+a) [l+O 923(ga? ^ ] 


13 9 

L 1 + 2.3*10 y 


dp 

1 25 

dt = 474.2. y. ( x-xra) 

Here 

xm An approximation to the begining of mach 

reflection region 

xo Approximation to the points in the double peak 
region where the peaks are equal, 
dt Time sepration between the first & second peak 
d Ratio of the first peak to second peak 
ta Approximation to the shock arrival time 
dp Approximation to the duration of the positive 

phase . 
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Clinical Effects of Ionizing Hadlailon Doses 


Apex ' c' 


Eftect'Phase 

0) 

DOSErRANGE 

0 to IDO remi 
— eubcllnice! 
range 

(2) 

100 to t,000 Rems— 
Therapeutic Range 

Over 1,000 Rems— 

Lethal Range 

100 to 

200 rema — 
cllnicat 
iurveltlance 
0) 

200 to 

600 rents— 
therapy 
effective 
( 4 ) 

600 to 

liDDO rema— 
therapy 
prorntslng 

(S) 

1,000 to 
5,000 rems— 
therapy 
paltiallye 
(6) 

Over 

6,000 rems— 
therapy 
palliative 
(7) 

lr>cidence of 
vomitir>g 

None 

100 rems 
infrequent 
200 rems 

common 

300 rerns 

100% 

100% 

100% 

100% 

tniliai phase 
Onset 

OuratKm 

— 

3 to 6 hr 
« 1 day 

VS to 6 hr 

1 to 2 days 

’4 to Vq hr 
c 2 days 

5 to 30 mir^ 
c 1 day 

Almost 

Unmediatel/* 

Latent phase 
Onset 

Dural K)n 

— 

« 1 day 
« 2 wk 

1 (0 2 days 

1 to 4 wk 

^ 2 days 

5 to 10 days 

« 1 day* 

0 to 7 days* 

Almost 

immediately** 

Final phase 
Onset 

Duration 

— 

10 to 14 days 
4wk 

1 to4 wk 

1 toe wk 

5 to todays 

1 to4 wk 

0 to ID days 
2 to 10 days 

Almost 

Immediately** 

Leading organ 

Hematopoietic tissue 

Gastrointestinal 

tract 

Central nervous 
system 

Characierlstfc 

stgns 

None below 
$0 rems 

Moderate 

ieukopenla 

Severe leukopenia; purpura, 
hemorrhage, tnfectian 

Epilation above 300 rems 

Diantiea; fever, 
disturbance of 
electrolyte 
baiarxDe 

Convulsions , tremor , 
afaxia, lethargy 

Crttical period 
postexposure 

— 

— 

1 toSwk 

Z (0 14 days 

1 to4£ hr 

Therapy 

Reassurance 

Reassurance, 

hematologic 

surveillance 

Blood trans- 
fusion, anti- 
biotics 

Conslderbone 
marrow trans- 
l^anlalion 

Maintenance of 

elecirolyfe 

balance 

Sedatives 

Prognosis 

Excetient 

Excellent 

Guarded 

Guarded 

Hopeless 

Convalescent 

perkxj 

Inddence of 
death 

None 

Norte 

Several wk 

None 

1 to 12 mo 

0to90% 

Long 

6010100% 

100% 

Death occurs 
within 


— 

2to12wk 

1 to6 wk 

2 to 14 days 

< 1 day to 2 days 

Cause o\ death 

— 


Hemorrhage, tnfectioo 

Ci/culalory 

collapse 

Respiratory failure, 
brain edema 


• At the highef doses within this range there may be no latent phase 

* InHiaJ phase merges Mto final phase, death osoally occumng from a lew hours to about 2 days, tfiis chronology Is possiWy 


Inlerrupled by a very short laleni phase. 


Apex D’ 

TAYLOR’S METHOD FOR CALCULATION OF BUILDUP FACTORS 

Taylor s forimuia for calculation of buildup factor is 


B = A. e + I 1 ~ A)e 


Where values ofA ,a, 13 for exposure buildup factor data 
for ordinary Concrete & soil of the composition given in Apex ‘E’ 
are given In the Tables below 

For Ordinary Concrete 

Photon Energy 1 A a 8 


0.058 


10.0 

For soil 


0 , 

5 

1 . 

,0 

2 . 

,0 

4 . 

,0 

6 . 

,0 

8 . 

.0 


D 056 
- 0.050 


0 084 


0.0754 

- 0.0471 

0.0653 

0.0224 

0.0538 

0.0599 

0 0388 

0,0690 

0.0357 

0.0675 

0.0421 

0.0920 
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